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Abstract. We study the equilibrium configuration of a nematic liquid crystal bounded by a 
rough surface. The wrinkling of the surface induces a partial melting in the degree of orientation. 
This softened region penetrates the bulk up to a length scale which turns out to coincide with the 
characteristic wave length of the corrugation. Within the boundary layer where the nematic degree 
of orientation decreases, the tilt angle steepens and gives rise to a nontrivial structure, that may be 
interpreted in terms of an effective weak anchoring potential. We determine how the effective surface 
extrapolation length is related to the microscopic anchoring parameters. We also analyze the crucial 
role played by the boundary conditions assumed on the degree of orientation. Quite different features 
emerge depending on whether they are Neumann- or Dirichlet-like. These features may be useful to 
ascertain experimentally how the degree of orientation interacts with an external boundary. 
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Nematic liquid crystals are fluid aggregates of elongated molecules. When the 
nematic rods interact with an external surface, the elastic strain energy induces them 
to align parallel to the unit normal, even if the surface is not perfectly flat pp. Recent 
experimental observations confirm that the surface alignment of the nematic director 
is completely determined by the roughness- induced surface anisotropy [2]. Further an- 
alytical calculations, performed within the classical Frank model with unequal elastic 
constants, have detected the bulk effects induced by a periodically-modeled external 
boundary 

A crucial effect, still related to surface roughness, escapes the framework of Frank 
theory, where the only order parameter is the director. Indeed, it is physically intu- 
itive that nematic molecules will disorder if we force them to follow a rapidly varying 
boundary condition. This surface melting was first experimentally detected in 0J|S]. 
Recent experimental observations have also measured a boundary-layer structure in 
the degree of orientation [SJ. The surface melting has been confirmed by approximated 
analytical solutions |7j, numerical calculations 8,9, and molecular Monte-Carlo sim- 
ulations |1U) . 

The combined effect of a rapidly-varying director anchoring and surface melting 
gives rise to an effective weak- anchoring effect that was first observed in The 
problem of relating the effective anchoring extrapolation length to the microscopic 
roughness parameters has been studied in several theoretical papers, all framed within 
the Frank theory ^1^3^]. This observation is of basic significance, since weak an- 
choring potentials have proven to influence deeply all nematic phenomena, beginning 
with Freedericksz transitions |151 1161 IT?] . Indeed, several theoretical studies have al- 
ready determined the influence on anchoring energies of the presence of permanent 
surface dipoles JHj or diluted surface potentials ^0 Ell- 
in this paper we analyze in analytical detail the boundary-layer structure induced 
by a rough surface which bounds a nematic liquid crystal. We frame within the 
Landau-de Gennes order-tensor theory, to be able to detect the effects on both the 
director and the degree of orientation. Our results confirm the surface melting already 
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obtained in but allow us to detect new phenomena. First, the nematic director 
steepens close to the boundary, so giving rise to an effective weak anchoring potential, 
that turns out to be deeply related to the surface- melting effect, and thus can be 
correctly handled only within the order-tensor theory. Furthermore, the boundary 
layers display a strong dependence on the type of boundary conditions imposed on 
the degree of orientation. Indeed, the orders of magnitude of all the expected effects 
depend on whether the boundary conditions are Neumann- or Dirichlet-like. We 
discuss how these effects may help in ascertaining in experiments how the mesoscopic 
parameter, which measures he degree of order, interacts with an external surface. 

The paper is organized as follows. In Section ^ we present the model, we set 
the geometry of a roughly-bounded sample, and derive the Euler-Lagrange partial 
differential equations that determine the equilibrium configurations. In Section El we 
perform the perturbative two-scales analysis that provides all the analytical details 
of the boundary-layer structure. In Section [21 we solve an effective problem, in which 
the rough surface is replaced by a weak- anchoring potential. The concluding Sec- 
tion m compares our outcomes with the effective results of Section |31 and draws the 
conclusions. 

1. Equilibrium configurations. The degree of order decrease has been recog- 
nized by many authors as a crucial effect of surface roughness 01111 • We thus describe 
nematic configurations in the framework of the Landau-de Gennes Q-tensor theory 
|2*T| . The order tensor is defined as the deviatoric part of the second- moment of the 
probability distribution of molecular orientations: 



where I denotes the identity tensor. Q is a second-order traceless symmetric tensor, 
with spQc [- §,§] E3- 

In order to keep computations simple, we adopt the one-constant approximation 
for the elastic part of the free energy functional 



where K is an average elastic constant. We stress, however, that it is straightforward 
to generalize all what follows to take into account unequal material elastic constants. 

The free-energy functional includes the Landau-de Gennes thermodynamic po- 
tential as well 



The material parameter A depends on the temperature, and in particular it becomes 
negative deep in the nematic phase. On the contrary, B, C can be assumed to be 
positive and temperature- independent. The potential (|1.3fl strongly favors uniaxial 
phases, in which at least two of the three eigenvalues of Q coincide. In fact, Q is 
expected to abandon uniaxiality mainly close to director singularities |22l 1231 We 
will not deal with any defect structure. Thus, though the uniaxiality constraint is not 
essential for our purposes, we follow the attitude of avoiding unnecessary complica- 
tions [251 I26| . and restrict our attention to uniaxial states 




(1.1) 



/ el [Q] = -K |VQ| 2 , 



(1.2) 



/LdG(Q) = AtrQ 2 - BtrQ 3 + CtrQ 4 . 



(1.3) 




(1.4) 
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The scalar s € [— \ , l] and the unit vector n are respectively the degree of orientation 
and the director. With the aid of (|1.4f) . the potentials (|1.2|) . l|1.3|) can be written as 

/oi[ S ,n] =if( S 2 |Vn| 2 + I|V S | 2 ) and / LdG (s) = |A s 2 - \B s 3 + \C S 4 . (1.5) 

When A < B 2 /(12C), the absolute minimum of the function /LdG(s) occurs at the 
preferred degree of orientation 



3B + \/9B 2 - mAC 
s P r := > . (1.6) 

In order to gain physical interpretation of the results, we also introduce the nematic 
coherence length £ and the dimcnsionless (positive) parameter u>, defined as 

e 2 :=^ and u 2 := - (s pI B - A A) . (1.7) 

The nematic coherence length compares the strength of the elastic and thermodynamic 
contributions to the free energy functional. We will show below that it characterizes 
the size of the domains where the degree of orientation may abandon its preferred 
value s pr . The number to depends on the depth of the potential well associated with 
s pr . Indeed, it is defined in such a way that /LdG( s pr) = uj2 /^ 2 - 

By using ljl.6 M l.7 p we write the bulk free-energy density /b := / e i + /LdG as 



/ b [s,n] _ , 2 , 1,_ , 2 , 1 ( a 4 



K = S 2 |Vn| 2 + -|V S | 2 + ^ ^ - -s* (2s pl - J- j + 2^(4, - co 2 ) j . (1.8) 

1.1. Modelling a rough surface. We aim at analyzing the effects that a rough 
boundary induces in a nematic liquid crystal. Once again, we try to keep our analysis 
as simple as possible, while still catching the essential features. We thus follow e.g. ^2] 
in modeling roughness by imposing a sinusoidally-perturbed homeotropic anchoring 
condition on a flat surface. The amplitude and the wave length characterizing the 
perturbation will be the crucial parameters in our results. 

We focus attention on a thin boundary layer, attached to the external surface. 
Consequently, we disregard the detailed structure of the bulk equilibrium configu- 
ration, that will only enter our results as asymptotic outer solution for the surface 
boundary layer. We introduce a Cartesian frame of reference {e x ,e y ,e z }, and assume 
that the nematic spreads in the whole half-space B = {z > 0}. We further simplify the 
geometry by assuming that n(r) = sm9(r)e x + cos#(r)e z and that the asymptotic 
bulk configuration depends only on z 

6(r) re 6 h (z) as z -> +oo . (1.9) 

In the presence of strong homeotropic anchoring on a flat surface, the boundary 
condition to be imposed on the director would be #( flat ) (x, y, 0) = . On the contrary, 
we will require 

e(x,y,Q) = Acos- . (1.10) 
V 

The boundary condition ijl.lUjl mimics the rugosity of the external surface by in- 
troducing two new parameters: the (dimcnsionless) roughness amplitude A and the 
roughness length r\. We remark that the oscillation rate increases as rj — > + , while 
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all roughness effects are expected to vanish in the limit A — > + . The requirements 
(|1.9|) . i|1.10|) imply that the free-energy minimizer will not exhibit any dependence 
on the transverse y-coordinate, so that we will henceforth restrict attention to the 
dependence on the coordinates (x,z). 

It is more complex to ascertain the correct type of boundary conditions which are 
to be imposed on the degree of orientation s. From the mathematical point of view, it 
would be natural to imitate the (Dirichlet) strong anchoring imposed on the director, 
and thus set s(x, y, 0) to be equal to some fixed boundary value s. Nevertheless, while 
it is well-established that we can induce an easy axis for the director on an external 
boundary, it is questionable whether we may fix the value of a mesoscopic parameter, 
that measures the degree of order in a distribution. From the physical point of view 
it would appear then more natural to impose (Neumann) free boundary conditions on 
the degree of orientation, leaving to the thermodynamic potential (|1.3fl the assignment 
of inducing the preferred value s pr in the bulk (z — > oo). To be safe, both possibilities 
(Dirichlet and Neumann) will be analyzed in Section [5] 

1.2. Euler-Lagrange equations. Once we consider that |Vn| 2 = |V#| 2 , it is 
easy to derive the Euler-Lagrange partial differential equations associated with the 
functional (|1.8(l . They read: 

s 2 A(9 + 2sVs- V6» = and As - 3 s I V0| 2 - 3 = , (1.11) 

where 

a(s) := s(s - s pr ) ( s - s pr + — ^ . (1.12) 

Since the boundary conditions (|1.10|) are x-periodic, with a period of 27r?7, we look for 
solutions of (|1.11|) in (the space of C 2 -functions, 27r?7-periodic in the x-direction) . 
To complete the differential system (|1.10|) . in H2.1l we will require 

0(x,O) = A cos- ( 0( s g() 

^ asz^oo (1.13) 



f^(x,0) = \s{x,z) 




. dz 

while in ij2.2l we will choose 

6(x,0) = Acos- 

V and "'^~> as z -> oo (1.14) 

s(x,0) 

2. Two-scales analysis. Before proceeding with the perturbation analysis of 
the differential equations, we state them in dimcnsionless form. It will turn out that 
the correct scaling is obtained by measuring lengths in 7y-units, so that we introduce 
the new dimensionless coordinates x = x/rj, z = z/f], and define the dimensionless 
parameter e — Equations (|1 . 1 1|1 become thus 

s 2 A6> + 2s Vs • V0 = and e 2 As - 3e 2 s |V6»| 2 - 3ct(s) = , (2.1) 



where both the gradient and the laplacian are now to be intended with respect to the 
scaled variables. The nematic coherence length is usually much smaller than all other 
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characteristic lengths. Consequently, we will look for uniformly asymptotic solutions 
to l|2.1|) . by treating e as a small parameter. In this limit, equation l|2.1|) 9 is singular, 
so that a regular asymptotic expansion would not provide a uniform approximation 
of the solution. Indeed, the small parameter e multiplies the highest derivative, so 
that we may expect the solution to have a steep behavior in a layer of thickness 
6 (to be determined), close to the boundary z = 0. We refer the reader to the 
books [23 1281 1291 ITTTT) for the details of the singular perturbation theory we will apply 
henceforth and, in particular, for the technique of the two-scales method which directly 
yields a uniform approximation of the solution. 

A standard dominant balance argument (that requires to introduce a stretched 
variable Z — z/6) allows to recognize that the boundary layer thickness is 6 = s. 
We then introduce the fast variable Z — z/e. The two-scales chain rule requires to 
replace d z by (d z + e~ 1 dz), and equations l|2.1|l take the form (when s ^ 0) 

s {e 2 9- x - x + e 2 9- z - z + 2e9 rzZ + 9 >zz ) + 2e 2 s x 6, s + 2(ea, g + ^z)(eB- z + <z ) = (2.2) 
e 2 s, ss + e 2 s 22 + 1es- zZ + s, zz - 3s [e 2 (0 s ) 2 + (ed- z + 6 tZ ) 2 } - 3a(s) = (2.3) 

where a comma denotes differentiation with respect to the indicated variable. In 
agreement with the two-scales method, 9 and s are now to be intended as 0[x, z, Z) 
and s{x,z,Z), that is functions of x,z and Z regarded as independent variables. It 
will be only at the very end of our calculations that we will recast the connection 
between z and Z: Z — z/e. We seek for solutions which may be given the formal 
expansions 

9(x, z, Z) = 6 {x, z, Z) + e0i(x, z, Z) + e 2 8 2 {x, z, Z) + 0{e 3 ) (2.4) 
s(x, z, Z) = sq(x, z, Z) + esi(x, z, Z) + e 2 S2(x, z, Z) + 0(e 3 ) . (2-5) 

If we insert ()2.4(I - (I2.5|I in l|2.2|l - H2.3|l . we obtain the following sequence of differential 
equations to O (1), O (e), and O (e 2 ) 



sq,zz - 3s (6» 0iZ ) 2 - 3ct(s ) = 

{40i,z) iZ + it (s 2 e ,z) z = -2 {s e 0>z ) - - 2s 0tZ e 0>z 

s\,zz - 6s 9o,z6i,z - 3si (<j'(s ) + {0 ,z) 2 ) = 6s o 0o,z#o,z _ 2s o,zz 
(i(4d2,z),z + j- 2 (sl9 , z ) tZ = -j-tfO^lz - ±{4e 0tg ), s - ±(s*9 , s ) 
- 2 (« o 0i,z) t8 ~ 2 (sxB QfZ ) ti ~ ^s hZ e ,s - 2a o ,z0i,g 

< &2,ZZ - 3S2 

+ 3s 



(2.6) 
(2.7) 



'JM + (0 ,z) 2 ] - Gs 9o, z d 2 ,z = ls 2 a"(s ) (2.8 
2 _l (a, \2 



,s + vi,z) +{V0,x) 
+ 60 ,z {s\9\, z + si# 0: 2 + s 6i^) - 2si : 2z - s ^ zz - s 0}XX . 

Analogous equations can be easily derived at any desired order. For any n > 1, the 
differential system obtained at O (e ra ) is linear in the unknowns n ,s n , and may be 
solved analytically. By virtue of the multiscale method, we find the correct dependence 
on z, Z by requiring that all s n , 9 n are uniformly bounded as e — > + for expanding 
intervals of the type < Z < Z*/e, for a suitable positive constant Z*. In most 
practical cases this requirement is equivalent to asking the removal of secular terms 
(i.e. terms that diverge as Z — » +oo). 
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2.1. Free surface degree of orientation. In terms of the scaled variables, the 
boundary conditions l|1.13fl require 

J*M>=*«»* and whenE>>v . (2 . 9) 

[s^(x,0)=0 ^s(x,z)ws pr 

We introduce m := 0' h (0), the derivative of the asymptotic solution at z = 0, since 
it will play an important role in the following discussion. The leading solutions in 
expansions (|2.4J) . (|2.5|1 are 

So(x,z) = St )I and 0q(x, z) — m z + Ae~ z ^ ri cos — . (2-10) 

V 

Higher order asymptotic solutions are gathered by means of laborious but straightfor- 
ward calculations. After recasting the solutions in terms of the dimensional variables 
x = rj x and z = rj z, we find 



8 ( X , z) = s pl - fm 2 - — cos - + ^ e" 2 ; " 

\ rj rj 2 

2 



+ HfE£L c -^/«( ^_ _^cos-)+ Ole 4 ) and (2.11) 

0(x,z)=mz + Ae- z /" C os- + ^ f 1 - e" 2 ^ 

rj uj * \ rj V 

-^-(e- z ^ ~e- 3z ^)cos-~^^ze~ z ^ cos-) +0(6*) . (2.12) 
2rj- \ J rj rj rj J 

The above expansions have been carried out up to the first nontrivial correction of 
the Oth-order approximation. Indeed, all calculations must be pushed to O (e 3 ) since 
an internal £-layer is necessary to satisfy the boundary condition l|1.13|l in z = 0. 
This layer is of O (e 3 ) because in the Neumann case the boundary condition l|1.13fl 
concerns the first derivative of s, instead of the degree of orientation itself. We remark 
that the solutions (|2.11|I - H2.12|I are coherently ordered for every fixed value of rj ^ 0. 
However, they are not uniformly ordered when rj € (0, 77] , namely we don't have a 
uniform solution if rj is allowed to become of order £ or, still worse, tend to zero. In 
other words, the above solutions remain valid as rj — * + if and only if £ = o (rj) . The 
main properties of the equilibrium configurations in the mathematically appealing 
but physically uncommon case in which 77 is of the order of, or even smaller than, £ 
will be presented elsewhere |.31| . 

2.1.1. Surface melting. We can highlight three different contributions in the 
degree of orientation l|2.11|l . First, we notice a uniform decrease in the degree of order, 
equal to — s pr m 2 £ 2 /w 2 . This disordering effect is triggered by the ^-derivative m, and 
was certainly to be expected. In fact, a glance to the free energy functional (|1.8fl 
suffices to show that a reduction in s decreases the free energy whenever the gradient 
of the director is not null. We then find two boundary layers. The former, of thickness 
rj and O (e 2 ), is a further reduction of the degree of orientation due to the boundary 
roughness, that induces a director variation in the cc-direction. An internal boundary 
layer, of thickness £ and order O (e 3 ) , is finally needed in order to cancel the normal 
derivative of s at the external surface. If we take into account all the contributions, 
the mean surface degree of orientation, defined as the x-average of s(x, 0), turns out 
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to be 



< s(x,0) > 3 



" pr 



m 2 e A a e 



2f2 



2 C 3 



9 9 



2A 2 £ 
V3 w 3 ?7 3 



(2.13) 



Figure 12.11 evidences the reported behaviour of the mean degree of orientation as a 
function of the distance from the surface. 



< s(x, z) > x 




v z/r) 

Fig. 2.1. Boundary layers in the mean degree of orientation < s(x,z) > x , when § = 0.25n, 
s P i = 0.8, ui = 0.6, m = O.l/n, and A = 1.5. The plot exhibits the presence of two boundary layers, 
the internal one being required by the free boundary condition applied on s. 



2.1.2. Effective surface angle. The tilt angle 8 exhibits a boundary-layer 
structure as well. Equation l|2.12|l shows that such layer is of O (e 2 ) and thickness rj. 
It gives rise to an interesting effective misalignment of the surface director. Indeed, 
if we allow z 3> -q in (|2.12() we find that 



(x,z) ~ 9h{z) 



2m£ 2 A 2 
•q io 2 



as z ^S> rj 



(2.14) 



The asymptotic approximation l|2.14|l shows that an experimental observation, per- 
formed sufficiently far from the external plate (with respect to the microscopic scale 
77) would detect an effective tilt angle $b, whose value at the plate is different from 
zero, since 



>(0) 



2m£ 2 A 2 

T]LU 2 



(2.15) 



Thus, a coarse observation of the nematic configuration measures a surface tilt angle 
slightly different from the homeotropic prescription 8 sur i = 0. Figure l2~2l evidences 
this effect. In the next section we will analyze in more detail the result (|2.15() . Then 
we will show how it matches the predictions of an effective weak anchoring potential. 
We remark that the tilt angle does not exhibit any further boundary layer at the 
smaller scale £. 

2.2. Fixed surface degree of orientation. The perturbative analysis of the 
differential equations (|1.11|) . with the Dirichlet boundary conditions (|1.14fl . would be 
unnecessarily entangled because of the non-linearity of the thermodynamic potential 




(|1.12(l . In fact, in this case only implicit solutions for So(x, z, Z) can be gathered. In 
order to pursue our analysis, and still catch the essential features of the solutions, 
we replace the function a in 11. 11( 1 by its linear approximation <Ji(s) = lo 2 (s — s pr ). 
This is tantamount to replacing the Landau-de Gennes potential in l|1.5[) by a tangent 
quadratic well, still centered in s pr . Such approximation is certainly well-justified 
deep in the nematic phase, when the isotropic state s = becomes unstable, and the 
second well of the Landau-de Gennes potential can be neglected. 

The asymptotic properties of the solutions in this case depend critically on the 
value s forced on the surface. If s ^ s pr , the boundary layer induced by the Dirich- 
let condition dominates over the roughness effect. Indeed, the leading asymptotic 
solutions are given by 



s(x, z) 



( Spr _ g) e -V^/C 

V3 (S pr 

3mA ( 1 - e~ z/r > 



s) i e -^ 

u> 



477 ' ' 



-2z/ V 



— m z 
2 



x A 2 f 
r\ 2i] V 



1 



WMcos — 



0(e 2 ) (2.16) 



9(x, z) = m z + 



Ae~ z/v cos - 

n 



V3l 



h 



-MO) 



e" 2 /" cos - 



where 



h(0 = log 



3 pr 



(s pr - s) e 



O (s 2 ) , 
(s pr - s) e~^^ 



(2.17) 



(2.18) 



determines the tilt angle variation within the boundary-layer. The bulk-asymptotic 
tilt angle is then given by 



0(x,z) 



>(*) = 



m£ 



1 °P r 
log-r- 

s 



5 pr 



as z~^> t) 



(2.19) 
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We remark that, when s =/= s pr , the leading contribution to 9b (0) is independent of A, 
and thus does not depend on the surface roughness. Furthermore, the effective surface 
tilt angle depends linearly on £, which makes it significantly larger than the prediction 
H2.15fl . derived with Neumann-like boundary conditions on s, which possesses an extra 
£/t7 (small) factor. Finally, we remark the fact that 9b (0) shares the sign of m if and 
only if s < s pr . We will return below on the physical origin and implications of this 
result. 

When the induced degree of orientation s does exactly coincide with s pr , all 
calculations simplify since ft.(£) = log s pr , and all first order correction in (|2.17|l vanish. 
We therefore push our perturbation analysis, and obtain 



s(x, z) 



9{x, z) = m z 
A 3 



A 2 
rf 



-2z/r, 



2mA 



- Z /V , 



rf 

e 

u> 



A 3 / 
2r? 2 V 



Ae - */" cos- 

-3z/r) 



n 

2mA 
V 

'2mA 2 



x 

cos — 
V 



+ 0(e 3 ) 



S(^(— ) 



x 

cos 



2 m 2 A 
V 



X 

COS — 



(2.20) 



(2.21) 



Equation l|2.21|) allows to compute the asymptotic tilt angle 9b, when s — s pr . In fact, 
once we drop all exponentially-decaying terms in (|2.21|) . we arrive at the interesting 
result that #b(z) does exactly coincide with (|2.14l) . that is with the expression we 
derived with a Neumann-like boundary condition on the degree of orientation. In 
fact, the complete expression (|2.21|) for the tilt angle 9(x, z) coincides with (|2.12|) up 
to O (e 3 )- Thus, any observation on the tilt angle is not able to distinguish among 
a free and a fixed boundary condition on the degree of orientation, as long as the 
imposed value s coincides with the preferred value s pr . This similarity between the 
Neumann and Dirichlet cases can be pursued further. Indeed, we can determine 
the O (e 2 ) -contributions in (|2.1(il) - (|2.17|l also when s ^ s pl . If we then use them to 
compute the O (e 2 ) -correction to the asymptotic tilt angle (|2.19|l . we arrive at the 
following expression, valid at O (e 2 ) for any value of s: 



8{x,z)^9b(z) = 
that yields 



m£ (, Spr . Spr-S 
-j=- log -4- + 

V3w \ s s 



2m£ 2 A 2 



7] OJ^ 



mz as z > -q, (2.22) 



9 b (0) = ^flog^ + V-i) + ^^ 



(2.23) 



The O (e 2 ) -contribution to the effective surface angle 9b (0) is thus fully a roughness 
effect, and does not depend at all on the type of boundary conditions imposed on s. 
On the other hand, equation H2.23fl confirms that the effective surface angle possesses 
also an O (e)-term when Dirichlet conditions are imposed on the degree of orientation, 
and s 7^ s pr . 

Figure I2"3l shows how the degree of orientation varies within the boundary layer, 
as s is fixed above, equal to, or below s pr . A double boundary-layer structure emerges. 
All plots exhibit a decrease of s in a region of characteristic size rj: this effect comes 
from the O (e 2 )-contribution. A similar surface melting was already presented and 
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discussed in Figure l2~Tl Close to the boundary, the O (l)-term proportional to (s — 
s pr ) e"^" 2 ^ settles the desired boundary value of s in a thin boundary layer of 
characteristic size £. 



< s(x, z) > 




1 1,2 3 

jj z/?7 

Fig. 2.3. Boundary layers in the mean degree of orientation < s(x,z) > x , when f = 0.25??, 
s P r = 0.8, = 0.6, m = O.ln, and A = 1.5, when Dirichlet-like boundary conditions are applied on 
the degree of orientation. The boundary degree of orientation s is respectively equal to 1 ( top ), s pT 
(middle), and 0.6 (bottom). 

3. Effective weak anchoring. Once the boundary layer effects fade away, the 
main macroscopic effect of a rough surface on the director orientation is to allow for 
an effective surface tilt angle #b(0), that apparently violates the homeotropic prescrip- 
tion 8(0) = (sec (I2.15|) and Ij2.23|l ). It appears then natural to check whether the 
same macroscopic effect may be modeled through a weak anchoring potential, acting 
on a smooth surface. In this section we pursue this similarity, and we derive a rela- 
tion connecting the microscopic roughness parameters with a macroscopic anchoring 
strength. 

To solve the weak-anchoring problem, we consider a nematic liquid crystal which 
still spreads in the half-space B = {z > 0}. To better compare our results with 
classical weak- anchoring models, we settle within Frank's director theory, and thus 
look for the equilibrium distribution that minimizes the free-energy functional 

F[n] :=K { \Vn\ 2 dv + W [ f w [n}da . (3.1) 
Jb Job 

The bulk free-energy density in the functional (|3.1|l can be derived from its order- 
tensor theory counterpart by setting s = 1 in 11.8fl . The anchoring potential / w is 
required to attain its minimum at the homeotropic anchoring n\ dB — e Z) while W is 
the anchoring strength. 

We look again for equilibrium distributions of the type n(z) = sin9(z)e x + 
cos 9(z) e z . Thus, the free energy functional H3.1J) per unit trasverse area can be 
written as 

f[6] :=K J 9' 2 (z)dz + Wf w (9(0)) , (3.2) 

where we assume /4(0) = and /^(0) > 0, in order to guarantee the homeotropic 
preference. The minimizers of (|3.2|l satisfy the trivial Euler-Lagrange equation 6" = 0, 
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and the boundary condition 

K9'(Q)-Wf^(9(0))=0. (3.3) 

When the anchoring strength W is large enough, the boundary condition 1)3.3(1 requires 
9(0) to be small. When this is the Taylor expansion in (|3.3I) supplies 

g(0) a Km = Am , (3.4) 

In 1)3. 4JI we have restored the notation m = 6>'(0) to better compare this estimate with 
our preceding results, and introduced the surface extrapolation length 

A := WfW) ' (3 ' 5) 

a quantity that compares the relative strengths of the elastic and anchoring potentials. 

The comparison between (I3.4|l and our results (|2.15l) - (|2.23|) relates the surface 
extrapolation length to the microscopic roughness parameters and/or the surface value 
of the degree of orientation. This analogy will be examined in the following section. 

4. Discussion. We have examined both the boundary layer structure and the 
bulk effects of a rough surface bounding a nematic liquid crystal. Our main results 
may be summarized as follows. 

• The roughness of the surface has been modeled by an oscillating anchoring 
condition, characterized by an oscillation amplitude A and a wave length rj. 
Figures IO and IO show that the rough boundary induces a partial melting 
in a neighborhood (of size rj) of the external boundary. When Neumann-like 
boundary conditions are imposed on the degree of orientation, equation 1(2.13(1 
quantifies the mean degree of order at the boundary. On the contrary, were s 
be forced to a prescribed value s on the surface, equations ((2. 161) and ((2.2U|) 
show that the boundary condition induces a thin boundary layer, determined 
by the nematic coherence length £. 

• Once the degree of orientation decreases, the spatial variations of the tilt 
angle become cheaper, and thus the 9 is keen to steepen close to the external 
boundary. Figure l2~2l illustrates this effect. As a consequence, the effective 
boundary tilt angle 6>b(0), extrapolated from the asymptotic outer solution 
9b(z) 7 becomes different from the null homeotropic prescription (see equations 
1(2.15(1 and ((2. 23(1 ). In the preceding section 13 we have shown that a similar 
effective anchoring breaking takes place when a weak anchoring potential is 
assumed on a smooth surface (see equation (13.5(1 for the characteristic surface 
extrapolation length). To further pursue this similarity we need to consider 
separately the different anchorings that may be applied on the degree of 
orientation. 

— When s is free to choose its boundary value, equation 1(2.15(1 shows that 
the surface extrapolation length is given by 

* £l +0 (£). ,4.1, 

Thus, the anchoring strength increases when either the roughness am- 
plitude A decreases (towards a smooth surface) or the roughness wave- 
length increases. An estimate of the order of magnitude of the effective 
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roughness wave-length can be obtained by assuming typical values for 
the quantities involved in l]4.ip . Indeed, if we assume A»{, Asl, and 
oj ss i we arrive at r\ « 10£, that models a roughness wave length in the 
hundredths of molecular lengths. 
— When the boundary conditions fix the value of the degree of orientation 
at the surface, equation 1)2.23)1 yields 

Equation 1)4. 2|l shows that the surface extrapolation length includes two 
quite different contributions. The former depends on the difference be- 
tween the boundary and the preferred values of the degree of orientation 
(s and s pr , respectively), while the latter depends on the surface rough- 
ness and indeed coincides with (|4.1J) . However, equation 1)4. 2JI may lose 
sense when s > s pr . Indeed, in this case A may become negative, so pro- 
viding an inverse weak anchoring effect. The physical origin of this odd 
result may be easily understood if we again resort to the s 2 |V#| 2 -term 
in the free-energy density. By virtue of that term, the tilt angle prefers 
to limit its spatial variations in regions of higher s. If we force in the 
surface a higher degree of orientation than the bulk value, the tilt angle 
will flatten close to the surface, thus exhibiting the opposite behaviour 
with respect to that shown in Figure l2~2l Equation l|4.2|l shows that this 
inverse effect may occur whenever 

l-±>^i + o(?). (4.3, 

If we again replace the estimates above for A,u>,T), we arrive at the 
result that a fixed degree of orientation is able to completely hide the 
roughness-induced effective weak anchoring whenever s exceeds s pr by 
the 10% of the preferred value s pr itself. 

Acknowledgements. P.B. thanks Georges E. Durand for useful discussions on 
the present topics. 
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